Topological Quantum Phase Transition in an = 2 Spin Chain 
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We construct a model Hamiltonian for S — 2 spin chain, where a variable parameter a is intro- 
duced. The edge spin is 5* = 1 for a = 0, and S = 3/2 for a = 1. Due to the topological distinction 
of the edge states, these two phases must be separated by one or several topological quantum phase 
transitions. We've studied the quantum phase transition by DMRG calculation, and proposed a 
phase diagram for this model. 
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I. INTRODUCTION 

Recently, investigations of topological phases and 
phase transitions has attracted great attention in con- 
densed matter physics The quantum Hall state [1] 
is the first example of a topological state of quantum 
matter, with a fully gap ground state in the bulk, and 
gapless excitations at the edge. The chiral edge state 
is a holographic mirror of bulk topology Q. In the 
recently discovered time reversal invariant topological 
insulators [3, 111, helical edge states are confined at the 
edge by the bulk energy gap, and states with opposite 
spins counter-propagate. In the case of the quantum 
spin Hall state realized in HgTe quantum wells, the topo- 
logically trivial and non-trivial states are separated by 
a topological quantum phase transition, tunable by the 
thickness of the quantum well. 

Quantum spin chain is another example where topolog- 
ical quantum phase transition is found. The low energy 
dynamics of ID large-spin Heisenberg antiferromagnet 
can be described as 0(3) nonlinear sigma modelQ. Half- 
integer-spin chains are generally gapless, while integer 
spin chains are gapped; they are described by the 0(3) 
nonlinear sigma model with and without the topologi- 
cal term. This distinction bears strong similarity to the 
topological insulators, which is distinct from the conven- 
tional insulators by the presence of a topological termQ. 

In this paper, we investigate a spin-2 chain model 
with two topologically distinct, transactionally invariant 
ground states. One model has edge state 5 = 1, while 
the other model has edge spin S = 3/2. The Berry's 
phase associated with the edge spins differ by tt. The 
square of the time reversal operator T gives = 1 for 
the first case, whereas it gives = — 1 in the second 
case. Due to this topological difference, the two ground 
states must be separated by one or several topological 
quantum phase transitions where the spin gap closes. 

This paper is constructed as follows. In the next sec- 
tion, we will review two exact solvable quantum spin 
model in one dimension. Some materials can be system- 
atically found in a brilliant work by Tu et.aljlj. After- 
wards, our new model Hamiltonian is presented accord- 



ing to the topological argument. Corresponding numer- 
ical results are shown in the third section, where phase 
diagram is discussed as well. Conclusions are drawn in 
the end. 



II. MODEL HAMILTONIAN 

Our starting point is an integer spin model introduced 
by Affleck, Kennedy, Lieb, and Tasaki, namely the AKLT 
modeljl, [l^. It is proved that the AKLT model has a 
unique infinite volume ground state, with an exponen- 
tial decay spin-spin correlation [loj. In agreement with 
the Haldane conjecture [ll|, the excitation gap of AKLT 
model is nonzero. The ground state of AKLT model can 
be written down exactly in terms of Schwinger bosons: 
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where S stands for site spin, and a\, b\ are creation op- 
erators of Schwinger bosons at the ith site. This ground 
state can be rephrased in a pictorial form. The spin-2 
AKLT ground state can be schematically shown in Fi^ 
where the circles stand for sites on the chain. Each spin-2 
can be decomposed into totally symmetric combinations 
of four spin- 1/2 states, and each state is represented by 
a solid dot in the figure. Two pairs of neighboring dots 
form singlet states, shown in red bonds. These bonds are 
usually called valence bonds, and in this sense, AKLT 
ground state is usually referred as Valence Bond State. 
It is proved rigorously that this ground state is unique 
under this periodic boundary condition. Due to the sym- 
metric intra-site coupling and anti-symmetric inter-site 
coupling, the parent Hamiltonian of this ground state is 
given by 
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where P3 and P4 are the projection operators onto the 
spin-3 and spin-4 subspaces respectively. The positive 
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FIG. 1: The sketch of spin-2 AKLT wavefunction. The circles 
represent sites, and each spin-2 is decomposed into four spin- 
1/2 solid dots. Red lines connecting the dots stand for singlet 
states. 

coefficients ensure the state in Eq.(IT]) to be the corre- 
sponding groundstate. 

In 1998, Scalapino, Zhang, and Hanke (SZH) intro- 
duced a SO (5) symmetric spin model 12] with an ex- 
act valence bond ground state. The original motiva- 
tion for the model is to illustrate the SO (5) theory of 
high Tc superconductivity, which unifies the antiferro- 
magnetic (AF) and the d-wave superconducting (SC) 
phases [Tsill^. However, it was soon found later that the 
S0(5) symmetry can also be interpreted as an enhanced 
spin rotational symmetry [l^. The SZH model contains 
five quantum states at each site, forming the vector repre- 
sentation of the SO (5) group. However, these five states 
can also be interpreted as the quantum states of the spin 
5" = 2 of the SO (3) spin chain. SZH presented an exact 
ground state wave function expressed as a matrix prod- 
uct state of the Dirac F matrices, and showed that the 
edge states of the SZH model are 4 fold degenerate at 
each edge. Interpreted as an S" = 2 chain language, the 
edge spin contains S = 3/2 spin quantum numbers. Fol- 
lowing this work, more valence bond states with higher 
symmetry groups have been constructed fsl. [l6l - [l9| . 

The basic idea of the SZH model is the following. The 
tensor product of two SO (5) spinor can be decomposed 
into S0(5) singlet, antisymmetric tensor, and symmetric 
traceless tensor, namely, 

5x5 = 1 + 10+14. (3) 

In analogy with the spin-1 AKLT model, the largest sub- 
space is projected out, leading to the desired SZH model: 

H = JY,Pu{xy), J>0. (4) 

(xy) 

Due to the Clifford algebra of the five F-matrics: F°F^ = 
25°'^ + 2iF°^, no symmetric traceless components are in- 
volved in the product of two F-matrices. As a conse- 
quence, 

1*"^^^) = Tr(F™iF™^ ...F™")|mim2...mAr), 

mi ^...mj^ 

(5) 

is the ground state of the above SZH model, where m 
is a vector label of the SO (5) group, which can also be 
interpreted as the rrii = —2,-1,0,1,2 quantum numbers 
of S* = 2 spin chain. 

Due to the relationships between the S0(5) and the 
S0(3) groups, there exists a natural deformation of the 
SZH model to an SU(2) spin-2 SZH model. The required 
map from S0(5) group onto SU(2) group with spin-2 is 



given by: 

10(S'O(5)) = 3(S'C'(3)) © 7(S'0(3)) (6) 
14(S'0(5)) = 5(50(3)) ©9(S'0(3)) (7) 

And therefore the SZH Hamiltonian deforms to; 

jjSZH ^ J2 J2P2(S„ S,) + J4P4(S„ S,), J2, J4 > 0. 

.(8) 

The ground state is unchanged up to an S0(5) rotation. 

To see if the SZH state is gapped or not, we can evalu- 
ate the ground state spin correlation function. The cor- 
relation of matrix product state can be easily derived by 
the transfer matrix technique, given by, 

(S'fS'^) = {TrG^)-^Tr[Z{S'')G'-^Z{S'')G^-''], (9) 

where /i = x,y,zi§l. Define g = X^m^™!"^)' then G = 
g''®g^Y.m r™ (8>r™, and = g''®S''g. As it's an 

isotropic magnet, the correlation functions are the same 
in any directions. After some detail calculation, we derive 
{SfS^) = (SfSy) = {SfS^} = -20 X 5-'- for integer 
r > 1. Therefore, the correlation length ^ of the SZH 
model equals to l/ln5. This finite correlation length 
indicates that the low-lying excitation in the SZH model 
is gapped, consistent with the Haldane conjecture. 

It is interesting to note that the correlation function 
of the SZH model is negative-definite, with a correla- 
tion length ^ = 1/lnh ~ 0.61. Consequently, the lat- 
tice constant is almost twice of the correlation length, 
and the spins at neighboring sites correlate extremely 
weakly. On the other hand, according to Arovas et.al's 
workf20'| , the correlation length for spin-2 AKLT is 1/ In 2 
which is roughly the lattice constant. Therefore, al- 
though AKLT model is also a strongly disordered anti- 
ferromagnet, the neighboring spins are closely correlated, 
leading to the conventional staggering correlation func- 
tion, say. Si -Sr cx (-l)'^2-^ 

Now we have two sets of models of 1-dimensional spin-2 
chain with exactly known Hamiltonians and ground state 
wavefunctions. The differences between the AKLT and 
SZH model are not only the analytic forms as they ap- 
pear, but also the topological distinctions. The same 
as topological insulator, the bulk topology is relevant 
to the edge state of an open chain. For the spin-2 
AKLT model, two solid dots at each edge remain free. 
Symmetrical combination of these two spin-1/2 dots re- 
sults an edge spin with S = 1. This boson- like edge 
state is consistent with large-N theory of SU(N) quantum 
antiferromagnets [2l| . However, the SZH model serves a 
complement of the large-N analysis. For an open chain, 
the SZH ground state is given by 

\^;i,j)= J2 (r™^r"^---r"").j|mim2...mjv). 

mi ,...m jv 

(10) 

It explicitly shows that at each edge, there are four de- 
grees of freedom since the matrix product state is four 
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dimensional. Therefore, the edge state of SZH model is 
spin-3/2, ie., fermion-like. That is completely different 
from the AKLT model. As the edge state is protected by 
topology, and is robust under perturbation, the AKLT 
and SZH models belong to different topological classes. 
It can be easily understood from the Berry phase's lan- 
guage. Berry phase $bp is the additional phase when 
the spin winds around, which relates to the expectation 
value of by exp(— i$Bp) = (T^), where T is the time 
reversal operator. It's well-known that = —1 for half 
integer spins such as S — 3/2, while = 1 for integer 
spins such as 5 = 1. As a consequence, the Berry phases 
of the two models under investigation differ by an angle 
of TT. It is this difference that makes topological distinc- 
tion of AKLT and SZH models. From another points of 
view, the edge spin determines the ground state degener- 
acy (GSD)[2j], and naturally serves as a topology index 
in analogy to fractional quantum Hall effect. These two 
models thus have different topology index due to different 
edge spins. 

Given the topological distinction of the two ground 
states, we construct a model Hamiltonian interpolating 
between the AKLT and SZH models: 



H{a) = (1 - a)H 



AKLT 
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SZH 



(11) 



Without loss of generality, we set J2 — K3 — 1, J4 — 
K4 = P in the following. As the edge state is robust un- 
less the gap closes, there must exist one or several topo- 
logical quantum phase transitions (TQPT) where the gap 
closes and reopens in the evolution of a from to 1 . This 
TQPT can be addressed by studying the behavior of en- 
ergy spectrum and correlation function at each a. 



III. NUMERICAL RESULTS 

The density matrix renormalization group (DMRG) 
method is emploved[23j in our study. For this purpose, 
it's helpful to rewrite the projection operators explicitly 
in terms of spin operators. Applying the identity 

2S ^ 

• S, = Y.k'^iJ + 1) - + l)]^jfe'), (12) 
j=o ^ 

one can easily get 
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— [162(S, . S,) - 7(S, . S,)2 - 10(S,; . S,)3 



-(s»-s,n + i 
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[90(S, • S,) + 63(S, • S,)2 + 14(S, • 

(15) 



2520'""' ' ^ 
+ (S, . S,f] 



For present study, we keep m = 600 — 1000 states in 
the DMRG block with more than 16 sweeps to get a 
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FIG. 2: (color online) Ground state energy per site Eo/N (a) 
and the absolute value of its corresponding second derivative 
\<fE/da''\ (b), as a function of a, got by DMRG with 600 
states and /3 = 1.0 at different system sizes. 



converged result, and the truncation error is less than 
10^^ in near the critical point, and much less than lO^^'' 
away from the critical point. We make use of the open 
boundary condition (OBC) and the total number of sites 
is iV = 600. To check the finite-size effect of the sys- 
tem, we have studied the ground state energy per site 
Eq/N and the absolute value of the corresponding sec- 
ond derivative \d'^E/da'^\ with respect to a, with system 
size TV = 10 - 1000 and /3 = 1, as shown in FigH Such 
the ground state energy Eq/N starts to converge, and 
the sharp peak of \d'^EQ/da'^\ appears at TV > 60. Be- 
sides the ground state energy, we have also calculated the 
correlation functions, and find that the finite-size effect 
can be neglected at TV > 600. Therefore, in the following 
calculation, N — 600 is adopted, for which the finite-size 
effect can be completely neglected. 

Consistent with the discussion above, the existence of 
TQPT is clearly supported by the numerical calculation, 
as shown in Fig. [2] The second derivative of the ground 
energy has a sharp peak around etc — 0.80, which in- 
dicates the second order quantum phase transition ad- 
dressed above. This is further confirmed by the real-space 
spin-spin correlation function, as shown in FiglS) At 
a = and 1, the correlation function is exactly the same 
with the analytical result. When a < 0.80, the correla- 
tion behaves similarly with the AKLT model, showing an 



4 



N=600 

i 


-B-DMRG 








(a)a=0.0 










J 


30 4( 



'(c)a=0.75 -0. 



ffiffi 



50 60 70 80 9Q 100 





lxlO"''r|rTTT- 

HI 




(e)a^0.8 


llilr^ 

SO 70 80 90 1( 



(g)a=0.85 .4x10- 



5 10 15 20 25 



20 40 60 




100 20 40 60 SO 100 



FIG. 3: (color online) The real space spin-spin correlation 
functions (Si ■ S, ) at different a, obtained by the DMRG 
with N = 600 and P = 1.0. 
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FIG. 4: (color online) Dimer order parameters of the model 
Hamiltonian with different choices of a and /3. (a)The exis- 
tence of dimerization is established by choosing different sys- 
tem sizes A'^. (b)Dimer order parameters for different p. This 
parameter disappears for the critical /3 — 0.27. 
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oscillating behavior with respect the lattice separation. 
However, when a > 0.80, the correlation behaves simi- 
larly with the SZH model, showing an negative-definite 
behavior. At the critical point, the correlation function 
undergoes an qualitative change from AKLT to SZH. 

To get the ground state phase diagram, we also cal- 
culate the peak position ac of {(PE/da^l as a function o.s 
of (3, as shown in FiglSj Above the red line, the system 
is topologically connected to SZH model, and therefore 
belongs to the same topological class as SZH phase. One 
feature of this phase is translational symmetric and has 
spin-3/2 excitations on the edge. Similarly, the regions 
where a is small belong to the same topological class as 
AKLT phase. Translational symmetry is also respected, 
but with spin-1 excitations on the edge. 

However, one thing we should keep in mind that even 
the model Hamiltonian is translational invariant, spon- 
taneous symmetry breaking is also possible. One well 
known example is the biliear-biquadratic model for spin- 
1 chain, whose Hamiltonian is also written in terms of lo- 
cal projections, and translational invariant. However, the 
emergence of the dimerized phase is addressed in previ- 
ous works 2^ 25 1 [i^. In order to quantitatively describe 
this issue, the following dimer order parameter 
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FIG. 5: (color online) Ground state phase diagram of the 
model Hamiltonian got by DMRG with A*' = 600. The loca- 
tion of second order phase transition for different f5 is indi- 
cated by the red line. The yellow region above red line belongs 
to SZH topological class, while the light blue region belongs 
to AKLT class. Possible dimerization is shown in the green 
region sandwiched between these two phases. 



Da — |(SiSi-|_i) — (Si+iSi+2)| 



(16) 



as a function of a is introduced, where i labels the center 
site in the spin chain so that one can minimize the possi- 
ble finite-size effect induced by open boundary condition. 
In the calculation, the total site number is set to be 
an even number to avoid potential ambiguity in the def- 
inition, in which case i is simply A^/2. The indications 
of dimerization is plotted in FiglH It is shown explicitly 
that dimerization appears with proper choices of parame- 
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ters a and (3. One can rule out possible finite size effect as 
the dimer order parameter doesn't scale with the system 
size N, shown in FiglD^a). When /3 decreases, the max- 
imum amplitude of dimer order parameter decreases as 
well, and the dimerization expands less and less regions of 
a monotonously. At the critical point j3 = 0.27, dimer or- 
der parameter vanishes for any a, and the corresponding 
dimerization phase shrinks as well. The system under- 
goes a quantum phase transition without spontaneous 
breaking of translational symmetry. During this phase 
transition, the energy and spin-spin correlation function 
have the same form as shown in Fig. ^ and Fig. ([3]). 
However it's worth mentioning that for Neel state of spin- 
2 chain, the corresponding dimer order parameter would 
be 8. Therefore our dimer order parameter is pretty small 
compared to strict antiferromagnets, and the dimeriza- 
tion phase cannot be confirmed definitely. One possibility 
is the gap between excitation state and ground state for 
the present model is too small to be distinguished nu- 
merically. As a result, certain dimerized excitation state 
enters into our results and lead to such finite but tiny 
dimer order parameter. A promising solution is to ap- 
ply the periodic boundary condition (PBC) here, so that 
one can rule out the possibility of dimerization acquired 
from open boundary condition. Therefore, we have also 
done some calculation by DMRG for system with PBC. 
For the system size N — 100 — 200 sites, we keep up to 
m = 3000 states with truncation error smaller than 10^'*. 
Finally, we find that both OBC and PBC systems give 
us consistent results. 

Observed the possible presence of dimerization, one 
can readily work out the phase diagram, see FiglS] The 
red curve stands for critical points of E / da^], while 
the two dashed dark curves are upper/lower bounds of 
dimerized phase of each /3 value. It shows the onset 
value for second order derivative coincides with the up- 



per bound of dimerized phase, so that they describe the 
same phase transition between dimerized phase and SZH 
phase. While on the AKLT side, the phase transition 
would be of higher order that is invisible in Fig 121 Despite 
the possible presence of dimerized phase, the system un- 
dergoes a quantum phase transition without spontaneous 
breaking of translational symmetry at ac when (3 = 0.27. 
As the symmetry is unchanged in this case, this phase 
transition is originated from the topology only. 

In conclusion, the topological distinction of the AKLT 
model and SZH model for the S — 2 spin chain is pre- 
sented in this work. A model Hamiltonian as an inter- 
polation between these two models is introduced. The 
quantum phase transition of the model is protected by 
topology, and established by DMRG calculation. The re- 
sults indicate the presence of a dimered phase between 
two topological phases of AKLT and SZH. This dimered 
phase shrinks at critical value of /3c — 0.27. One would 
like to realize this topological phase transition in real ma- 
terials. On the AKLT side, we have a well defined exam- 
ple already. S = 1 edge spin is observed in 5* = 2 chain 
of CsCri_2,Mga;Cl3[23. However, examples on SZH side 
are still missing. Probably one can employ cold atom 
techniques to realize TQPT proposed in this paper. 
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